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Abstract. To investigate the presence of small scale structure in the spatial distribution of H2 molecules we have 
undertaken repeated FUSE UV observations of the runaway 09.5V star, HD 34078. In this paper we present 
five spectra obtained between January 2000 and October 2002. These observations reveal an unexpectedly large 
amount of highly excited H2. Column densities for H2 levels from (v = 0, J = 0) up to (v = 0, J = 11) and for 
several v = 1 and v = 2 levels are determined. 

These results are interpreted in the frame of a model involving essentially two components: i) a foreground cloud 
(unaffected by HD 34078) responsible for the H2 (J = 0, 1), CI, CH, CH+ and CO absorptions; ii) a dense 
layer of gas (n ~ lO** cm~^) close to the O star and strongly illuminated by its UV fiux which accounts for the 
presence of highly excited H2. Our model successfully reproduces the H2 excitation, the CI fine-structure level 
populations as well as the CH, CH^ and CO column densities. 

We also examine the time variability of H2 absorption lines tracing each of these two components. From the 
stability of the J = 0, 1 and 2 damped H2 profiles we infer a 3a upper limit on column density variations 
AN(H2)/N(H2) of 5% over scales ranging from 5 to 50 AU. This result clearly rules out any pronounced 
ubiquitous small scale density structure of the kind apparently seen in HI. The lines from highly excited gas are 
also quite stable (equivalent to AN/N < 30 %) indicating i) that the ambient gas through which HD 34078 is 
moving is relatively uniform and ii) that the gas flow along the shocked layer is not subject to marked instabilities. 
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1. Introduction 

Studies of the molecular phase of the interstellar medium 
are severely hampered by the fact that the most abun- 
dant species, H2, cannot be observed from the ground. 
"Tracers", like CH, CO and its isotopomers are often used 
as surrogates; however, because of selective photodissoci- 
ation near the illuminated boundaries of clouds, deple- 
tion onto grains in their depths or, simply, time evolution 
under the coupled effect of internal dynamics and chcm- 
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* Based on observations performed by the FUSE mission and 
at the CFHT telescope 



istry, no tight relation is expected between the local abun- 
dance of H2 and those of minor species. In other words, 
the widely used notion of " standard" ratios can hardly be 
justified, even if, when integrated along the line of sight, 
the amounts of various species appear relatively well cor- 
related. 

The above limitation is especially severe when we 
deal with small scale structure studies. Within translu- 
cent molecular clouds, H2CO, OH and HCO+ molecules 
have been found to display AU-scale structure (Moore & 
Marscher 119951 Liszt & Lucas ODOOl) . In contrast, the dis- 
tribution of dust grains is quite smooth at small scales, 
both in translucent and dense material (Thoraval et al. 
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UnHnmSZl IIM3 Lada et al.HnSI Alves et aL l^IMIl . If 

H2 were to display marked AU-scale structure, as minor 
molecular species apparently do, then the smooth behav- 
ior of dust grains could be understood as a result of their 
large inertia (Thoraval et al. 11999(1 . However, it could also 
be that the H2 distribution is smooth at small scales with 
no or little density structure (like dust grains) and that 
"chemical" inhomogeneities are present. Clearly, direct in- 
vestigations of the distribution of H2 itself are needed to 
clarify the above issues and determine whether the struc- 
ture seen for minor species is due to fluctuations of their 
local abundance or rather reflects real density structure. 

We have thus taken advantage of the unique oppor- 
tunity offered by the successful launch of FUSE and pro- 
posed a GI program which aims at probing the distribu- 
tion of H2 over spatial scales ranging from about 1 to 
100 AU. A bright 09.5V star, HD 34078 (AE Aur) has 
been selected for repeated observations and, thanks to its 
large proper motion (Moffat et al. 11998(1 . information on 
the structure of the foreground medium can be obtained 
within a relatively small time interval. Indeed, at a dis- 
tance of 530 pc the implied transverse velocity is 103 km 
s"^ or 22 AU/yr. Its color excess, E(B-V) = 0.52, is such 
that a significant fraction of the foreground gas should be 
in molecular form. It is also of interest that the velocity 
structure of the molecular gas along this line of sight is 
quite simple, with only one major component that has b 
« 3 km s-i (Allen IT^ RoUinde et al. lMTSji . In the fol- 
lowing, we shall adopt a distance of 380 pc for this molec- 
ular material (Brown et al. 11995(1 . implying a drift of the 
hue of sight through the cloud at a velocity of 17 AU/yr. 
Note finally that HD 34078 has been the subject of sev- 
eral detailed studies in the past decade that have brought 
results on a number of molecular species including CH, 
CH+, CN, C2 (Federman et al. IBMl Herbig[B99j. The 
latter can be helpful for modelling and the past measure- 
ments provides information on structure over larger scales 
by comparison with present-day observations. In partic- 
ular, a series of observations of the blue CH and CH"'' 
lines has been performed since 1999 in parallel with this 
FUSE program (RoUinde et al. 12003(1 . Results indicate a 
long-term increase of N(CH) over the past 10 years while 
N(CH+) remained stable or decreased. Correlated short- 
term (~ 1 yr) CH and CH+ fluctuations might also be 
present. 

The first FUSE spectrum obtained revealed large 
amounts of high excitation H2. This is an unexpected find- 
ing in comparison with observations of most other lines 
of sight (see however Meyer et al. 120011 for HD 37903 and 
Federman et al. ll995(l . As HD 34078 is a high velocity star, 
it is likely to be far away from its parent cloud and should 
hence not be associated with molecular material. This re- 
sult has important implications for the identification of 
physical processes responsible for the unusual excitation 
but could be a difficulty for our structure study. Indeed, 
the presence of highly excited H2 raises the possibility of a 
significant contribution to the observed absorptions from 
gas associated with the star itself, which could be subject 



Table 1. Properties of the line of sight. Right side of the 
chart: extinction curve coefficients. 
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(1) Simbad database 

(2) Crawford tl975> 

(3) Patriarchi et al. II200H 

(4) Fitzpatrick and Massa I1990II 



to variability of a different nature than the one due to 
structure in the foreground gas. The first objective of this 
paper therefore becomes to establish a reasonable scenario 
for the origin of the observed absorption, based on detailed 
modeUing of the HD 34078 FUSE spectra. Secondly, we 
present a study of the variability of the H2 lines based on 
five spectra spanning a period of 2.7 years. 

In Sect. 2, we describe the FUSE observations and data 
analysis. Complementary OH observations performed at 
the CFHT are presented in Sect. 3. Results on the gas 
properties are discussed in Sect. 4 with their implications, 
as inferred from modelling of the physical conditions in the 
absorbing material. We next investigate H2 line variations 
(Sect. 5). Preliminary analyses of the first three FUSE 
spectra have been presented by Le Petit et al. ((20001 1^01(1 
and Boisse et al. ((2001(1 . 

2. The FUSE data 

2.1. Observations and data reduction 

The first FUSE spectrum was obtained on 2000 January 
23; it consisted of eleven exposures with a total integration 
time of 7145 s . Four other spectra were acquired in com- 
parable conditions on 2000 October 30, 2001 February 16, 
2001 October 15 and 2002 October 10 (cycle 1 and 2 pro- 
grams). The data were taken in histogram (HIST) mode. 
In this paper we rely mostly on the first spectrum for the 
study of the intervening gas properties, but occasionally 
we use the other four, either to secure detections of faint 
lines or to improve the accuracy of column density esti- 
mates. The search for variable H2 absorption has been 
performed mainly by comparing the 1st and 5th spectra, 
because i) they turn out to be of very similar "quality" 
(resolution and S/N), ii) they sample the largest time in- 
terval (i.e. largest spatial scales, over which changes might 
be easier to detect in the event of a regular variation) . For 
absorption lines that are suspected to vary in time, mea- 
surements are performed on the first spectrum. 

Standard processing using the FUSE pipeline software 
was applied to perform the wavelength and flux calibration 
(for details on the FUSE mission, see Moos et al. 120001 and 
Sahnow et al. i2000,l . The analysis presented here is based 
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mostly on the LiFlA and LiFlB spectra, which display 
the highest S/N ratio. 

The data were reduced in a standard way us 
ing MIDAS, the ESO data reduction software package 
Comparison of individual exposures occasionally revealec 
variations in absorption line positions from one exposur( 
to another. Nevertheless, these remain quite small and th( 
average spectrum obtained after applying the appropriate 
shifts does not differ significantly from the raw average 
The original pixel corresponds to 0.005 A but, given th( 
resolution provided by FUSE (about 20 km s~^) the spec 
tra were binned to 0.015 A. In the co-added spectrum 
the S/N ratio reaches 30 per 15 mA pixel which translate: 
into a 3 (T limit on the equivalent width (W) as low as ' 
mA. 

2.2. Line identification and column density 
determinations 

Line identification has been performed using the list pro- 
vided by Morton H2000(l for atoms and ions. For H2 tran- 
sitions we use the data from Abgrall et al. (1993ab). 
Originally, only lines from v = with J = 0, 1, ....9 levels 
were considered. As it turned out that lines from higher 
excited levels were present in our spectra, the table was 
complemented by adding transitions from v = 0, J = 10, 
11, 12, 13; V = 1, J = 0, 9 and V = 2, 3, 4, J = 0, 1, 2. 

To perform the line identification, we proceed in two 
steps. We first measure the position of all prominent ab- 
sorption lines in the spectrum and assign one (or several) 
transitions to each of them. Next, we search in a system- 
atic way for all strong transitions expected in our range 
and examine whether all detections/non-detections and 
W values for any given species (or particular level of a 
species) are mutually consistent. This step is a crucial one 
here since, given the high number density of absorption 
lines in the spectrum, blends are numerous and the con- 
tinuum level is often ill-defined (we estimate that in the 
1040 - 1060 A interval for instance, the fraction of the spec- 
trum free of absorption is no larger than 25%; this fraction 
becomes close to zero above 1100 A). As a result of this 
crowding, some lines may appear either too strong in the 
event of an unresolved blend or too faint if adjacent ab- 
sorption leads to an underestimate of the continuum level. 
When an unresolved blend or ill-defined continuum is sus- 
pected, the corresponding feature is not used for column 
density determination. A portion of the LiFlA spectrum 
is shown in Fig. ^ where several high excitation II2 lines 
are identified. 

Except for damped II2 lines, most features are not re- 
solved (recall that b ~ 3 km s~^ for molecular gas towards 
HD 34078). Further, the line spread function is not accu- 
rately known and may vary across the spectrum; we then 
often rely on curve of growth analysis to derive column 
densities. However, line fitting is required to deal with par- 
tially resolved blends and damped II2 lines. To this pur- 
pose, we use the routine OWENS developed by Martin 




1074 1076 1078 1080 

Wavelength (A) 



Fig. 2. Portion of the October 2000 LiFlA FUSE spec- 
trum comprising the strong II2 (2-0) Lyman band around 
1078 A (red line). The fit obtained when including the 
nine II2 lines indicated (v, J values for the lower level are 
given) is superimposed (thin line); thick tick marks indi- 
cate the four strong damped H2 lines from J = 0, 1, 2. 
The absorption near 1076 A is from CO. 



Lemoine and the LYMAN context in MIDAS. In a first 
step, the spectrum is shifted to zero velocity using narrow 
features located blueward and redward of the absorption 
considered. When fitting damped II2 profiles, the contin- 
uum is interpolated using low order polynomials. Good 
fits are easily obtained as shown in Fig. 12 For low J val- 
ues, several II2 systems are available that can be used to 
get independent column density estimates (hereafter N), 
which helps to further reduce the uncertainties. 

Damped lines from low energy II2 levels provide a mea- 
surement of N which is independent of the velocity distri- 
bution of the gas; for a b parameter of about 3 km s^^ (as 
measured for CH), we find that this is true up to J = 3. 
The situation is less comfortable for lines of intermediate 
optical depth. Indeed, difficulties in defining the contin- 
uum or due to blends are such that for most II2 levels we 
cannot build a curve of growth appropriate for deriving 
both the b parameter and N. We have to assume some a 
priori b value. To investigate the possibility of a variation 
in b with excitation energy, we searched for highly excited 
levels (J > 5) for which our data can be used to constrain 
the b value. Thus, for the (v = 0, J = 9) level, we get re- 
liable measurements for 8 transitions with log(A f) values 
ranging from 0.75 to 1.57 (cf. Fig.|2|). We also include three 
measurements from the (v = 1, J = 3) level which is close 
to the previous one in the energy scale. As can be seen in 
Fig. 01 both levels give consistent constraints on the curve 
of growth. From these, we estimate that a b value of about 
4 km s^^ is consistent with the data and find that values 
below 2 and above 6 km s^^ are excluded. We therefore 
find no evidence that the velocity distribution of the high 
excitation H2 differs from that of CH, which presumably 
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Fig. 1. Portion of the January 2000 LiFlA spectrum in units of 10^^^ erg s^"'^ cm~^A^^ near strong H2 (J = 0, 1, 
2) lines from the (4-0) Lyman band, v and J values for the lower level of each H2 transition involved are given. The 
features indicated illustrate the variety of excited H2 levels from which absortion is detected. The identification of the 
feature noted "4, 2 ?" is uncertain (see text in Sect. 2.3). 



reflects that of the bulk of H2. Note that Meyer et al. iPMl 
reach a similar conclusion for highly excited H2 detected 
in front of HD 37903 (the b value being 1.8 km s^^ in this 
case). In practice we shall adopt b = 4 km s^^ for J > 3. 

Optically thin H2 lines can also be of interest for deriv- 
ing N without any assumption on the velocity distribution, 
and for providing reliable constraints on the upper end of 
the excitation diagram. For b = 4 km s~^we estimate that 
a line around 1000 A and with an equivalent width of less 
than 8 mA can be considered to be optically thin (the N 
values derived either from the exact N(W) relation or the 
linear approximation differ by less than 10%). This is sig- 
nificantly above our detection limit and we do detect H2 
features of this type from the v = 0, J = 10 and J = 11 
or some v = 1 levels. Non-detections in regions with good 
S/N can also be used to derive upper bounds on N (e.g. 
for H2 v=0, J = 12 and 13). 

Table IT^ summarizes the results obained for N(H2, v, 
J) and gives the best values together with 3 cr lower and 
upper bounds. The latter have been estimated in a con- 
servative way. When different lines from the same level 



provide discordant results for N, we adopt the best com- 
promise by weighting various estimates according to the 
presence of suspected blends and quality of the data. The 
number of lines effectively used in the determination of N 
IS given m Table 10 as M/. For lines which are neither in 
the damped or optically thin regime, an uncertainty on b 
of ± 1 km s~^ was considered. 



2.3. Comments on specific H2 levels or species 

In the following, we give some details concerning the col- 
umn density determination for transitions: 

• II2 J = 0, 1, 2 : For the absorption from these levels, 
the continuum level and shape across the broad profiles is 
difficult to define. To get an estimate of the unabsorbed 
spectrum, we have artificially removed all narrow features 
other than those from H2 J = 0, 1, 2. We then simultane- 
ously fit the latter using OWENS; a polynomial form of 
order 3 has been adopted for the continuum. Fortunately, 
N(J = 0, 1, 2) is well constrained by the profile near the 
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Fig. 3. Curve of growth for H2 transitions from the v = 0, 
J = 9 and the v = 1, J = 3 levels (filled and empty triangles 
respectively). The horizontal scale is appropriate for b = 
4 km s^^; the slopes of the curve of growth corresponding 
to b = 2 and b = 6 km s~^ are indicated (dotted and 
dashed line, respectively). A few representative 1 a error 
bars are given. 



bottom of the line and weakly depends on the exact shape 
of the adopted continuum. 

• H2 J = 3, 4, 9 : for these levels, the number of 
lines available range from 1 (J=7) to 10 (J=9). 

• H2V = 0,J = 11: This is the highest v = transition 
detected. Two features are seen at 1045.22 and 1045.35 A 
with W = 10 ± 2 mA and W = 7 ± 2 mA. 

• H2 V = 0, J = 12 and 13 : Lines at 1057.62 and 
1072.73 A from J — 12 are expected in good "windows"; 
from their absence we get the upper limit N < 2.4 10^^ 
cni^^. The strongest transition in our range from J = 13 
is expected at 1040.70 A (f = 0.0237) and is marginally 
present (2 a detection). We infer the upper limit: N < 2.7 
10^3 cm-2. 

• H2 V = 1 : several features are seen from each of 
the J = to 5 levels. Therefore, their identification is 
secure. In contrast, no feature from the (v = 1, J = 6, 7, 
8) levels could be measured due to blending with strong 
absorptions. 

• H2 V 2, 3, 4 : A blend of the 1046.62 and 1046.66 A 
lines (from v = 2 , J = 1 and respectively) is marginally 
detected with W = 12 mA in the 1st, 4th and 5th spectra, 
which have the best resolution and S/N in this region. 
Unfortunately, we do not find other lines suitable for the 
determination of N(v = 2 , J = 0) and N(v = 2 , J = 1) 
separately. 

Regarding v = 3 levels, several lines are suspected to 
be present but have a low significance level. For instance, 
a line from (v = 3, J = 1) is seen in the 1st and 4th spectra 
at 1116.405 A with W = 3 ± 1.5 mA. The corresponding 
column density is N = 6 lO^^cm"^. However, most of these 
lines are expected at A > 1100 A where the crowding is 



quite severe, and we prefer to retain only upper limits for 
these levels. 

In an earlier report (Boisse et al. 1200 l|l . we assigned 
the 1046.45 A feature to the v = 4, J = 2 level (hue 5 in 
Fig. 1). This line appears clearly only in the first spectrum 
and given i) the small f value (6 10~^), ii) absence of other 
detections from the same level and iii) absence of time 
variations of clearly detected highly excited H2 features, 
this identification appears questionable. Only spectra from 
the HST/STIS would allow us to unambiguously quantify 
the amount of H2 in levels with excitation energies above 
10^ K. 

• HI : In our data, the Ly/3 line is present although 
heavily blended with H2 J = 0, 1, 2, 5 absorptions near 
1025 A (Ly7 and Ly(5 can also be seen but the S/N ratio 
degrades rapidly at lower wavelengths and blends are nu- 
merous so we use only Ly/3 to constrain N(HI)). Thanks to 
the large width of this Ly/3-H2 absorption, airglow emis- 
sion does not affect its profile (emission occurs only in the 
central part of the trough, where the flux is zero; Fig. lll|l . 
To fit the spectrum in this region, we adopt the values for 
N(H2, v=0, J) inferred from other H2 systems and vary 
N(HI). The shape of the red edge of the broad 1025 A 
trough (around 1027 A) is primarily due to Ly/3 (see Fig. 
[TT|l from which we estimate N(HI) — 2.7 x 10^^ cm~^, a 
value significantly larger than that derived by Mc Lachlan 
& NandyfTSS¥(N(HI) = 1.75 x lO^i cm'^) from lUE Lya 
data taken in 1979. We come back to N(HI) variations in 
Sect 5.4. 

• HD, J = : Several lines are detected (near 1007, 
1011, 1021, 1031, 1042, 1054, 1066 A) but unfortunately, 
most of them are either blended, located in regions of poor 
S/N or have an ill-defined continuum. For instance, the 
strong 1021.456 A line is expected to be blended with 
an SIX feature with an unknown oscillator strength. The 
1031 A feature is contaminated by an H2 (v = 1, J = 3) 
line for which N is reasonably well determined; OVIA1031 
might also contribute, but this cannot be assessed since 
the OVIA1037 line coincides with a damped H2 line. The 
1054 A line appears much too faint with respect to the oth- 
ers (W w 15 mA); probably, this is due to an ill-defined 
continuum caused by extra absorption adjacent to this 
line. In contrast, the 1066.27 A feature is too strong with 
respect to the 1031.91 A one. It is thus difficult to infer a 
value for N(HD, J = 0) consistent with all lines detected. 
Moreover, N values are strongly dependent on the veloc- 
ity distribution, which might differ from that of CH. We 
therefore cannot do better than give an order of magni- 
tude estimate: N(HD, J = 0) w lO^^-^^ cm~2. 

• HD J = 1 : Lines from excited HD have already been 
detected in the Copernicus spectrum of ^ Oph (Wright & 
Morton 1979)1 . The strongest feature mentioned by these 
authors, at 1021.916 A, is apparently present in our five 
spectra, although blended with other weak lines. We esti- 
mate ly = 8 ± 2 mA. Other lines of comparable strength 
are expected from HD J = 1; at the position of the 
1054.722 A transition, a strong line with W « 35 mA is 
present but it is much too strong in comparison to the pre- 
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vious feature and is likely of stellar origin. Finally, around 
1043.288 A (on the blue wing of a H2 J = 3 line) the spec- 
tra display large variations which are of instrumental or 
stellar origin and the search is inconclusive. The tentative 
detection of the 1021.916 A hne implies N(HD, J ^ 1) 
5.4 ±1.3 X 10^3 cm-2. 

• CI : Numerous lines from CI and its fine structure 
excited levels are detected which can be used to mea- 
sure N(CI), N(Cr) and N(Cr*) and get constraints on 
the density, n. This key parameter has been estimated by 
Federman et al. (1994) from the analysis of 2 excitation. 

Several combinations of lines from three multiplets 
were considered between 1129 and 1160 A; spectral re- 
gions containing lines with unknown f values were ex- 
cluded. Various attempts were made to fit line profiles us- 
ing OWENS; we successively considered either unblended 
or unsaturated features alone or all of them together. The 
instrumental LSF was allowed to vary slightly across the 
spectrum and b values between 3 and 6 km s~^ were con- 
sidered. The results are the following: 

- N(CI) = 9.4 (+ 7.0, - 5.8) 10^^ cm'^, 

- N(CI*) = 5.8 (-f 0.0, - 4.2) IQis cm^^, 

- N(Cr*) = 2.2 (+ 1.8, - 1.1) 10^5 cm-2. 

These estimates and uncertainties represent the best 
compromise we find given our limited knowledge of the 
LSF and velocity distribution. 



3. The CFHT data 

We observed OH absorption lines around 3080 A from 
the A-X (0-0) transition band towards HD 34078 with 
the Canadian-France-Hawaiian telescope. The observa- 
tions were performed in 2001, October 4-7. We used the 
GECKO spectrometer, in the UV f/4 Coude train at a 
nominal resolving power of R = 120 000. The CFHT EEVl 
CCD was used as the detector. The spectrum was cross 
dispersed by a grism and OH lines appear in the 18th or- 
der. An integration time of 14 hours provided a S/N ratio 
of 36 per 35 mA pixel. Data were reduced in a standard 
way using the IRAF package. The three OH lines expected 
in the range covered are clearly detected. An extract is pre- 
sented in Fig. 01 and equivalent widths are given in Table 
El The inferred OH column density towards HD 34078 is 
N(OH) = 3.5 ± 0.5 xlOi3 cm-^. 



4. Properties of the gas along the line of sight 

In this section, we study physical properties such as exci- 
tation, temperature and kinematics. Small scale structure 
is considered in the next section. We first discuss parame- 
ter values as inferred from the analysis of absorption lines 
measured in the first spectrum and then present a model 
attempting to account for all known observational con- 
straints. 



1.05 



0.95 



HD 34078 
CFHT GECKO 




3079 3080 3081 

Fig. 4. Portion of the normalised CFHT spectrum com- 
prising two OH absorption lines (indicated by tick marks). 
The broad absorption around 3079.9 A is due to a detector 
defect. High intensity values on two "hot pixels" (noted 
*) have been cut. 
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Fig. 5. Observed H2 excitation diagram (filled squares 
with 3cr error bars shown; above 11000 K, observed values 
are upper limits)) and predictions of the model. The en- 
ergies of the lowest vibrationally excited v=l, 2 H2 levels 
are indicated. 



4.1. Observational results and comparison to other 
lines of sight 

In Fig. Owe give the excitation diagram corresponding to 
the estimates or upper limits listed in Table 

As observed for other lines of sight, there is a clear 
flattening as excitation energy increases. The J = and 1 
values provide a direct measurement of the gas tempera- 
ture, T = 77 K (71 - 83 K). The slope inferred from the 
J = 3 to 7 values correspond to Texc — 360 K, whereas 
the detections up to energies of 10^ K are roughly fitted 
by Texc — 1200 K. No J parity effect is noticeable. 

In order to assess whether the H2 excitation towards 
HD 34078 differs from the one determined on other lines 



p. Boisse et al.: A far UV study of interstellar gas towards HD34078 



7 





10° 




10-^ 




10-2 




10-^ 




z 

ti 


10-^ 




10-5 




10-^ 








10-^ 




10-^ 




10-^ 



HD 34078 
HD 37903 

HD 102065 
HD 108927 
?Oph 
^ Per 



1000 



2000 
Energy (K) 



3000 



4000 



Fig. 6. Comparison of the H2 excitation diagrams ob- 
served for the hnes of sight towards HD 34078, HD 37903, 
HD 96675, HD 102065, HD 108927, C Oph and C Per. 
Column densities for the J = to 6 levels have been nor- 
malised to the total H2 column density for each line of 
sight. For HD 102065, the N(J) values correspond to a b 
value of 3 km s~^. 



& Nandv IT 984(1 . The corresponding column densities are 
summarized in Tablc0| It appears that, on one hand, col- 
umn densities of molecular species as well as the H2 tem- 
perature, Toi, are qualitatively similar to those found on 
other lines of sight and, on the other hand, the H2 exci- 
tation above J = 2 is unusually high compared to that in 
standard diffuse ISM. This suggests the presence of two 
distinct components along the line of sight: i) "quiescent" 
gas within a translucent cloud, responsible for the absorp- 
tion due to H2 J = and 1 and to other molecular species 
and ii) gas located close to HD 34078 dominating absorp- 
tion from H2 at J > 2. 



The HI/H2 ratio only provides a lower bound for com- 
parison with model predictions since diffuse HI not asso- 
ciated to molecular gas is likely to be present along the 
line of sight. Adopting the N(HI) value derived from our 
FUSE observations, the resulting total hydrogen column 
density along the hne of sight, N(H) = N(HI) + 2 N(H2), 
is 3.3 • 10^1 cm-^. Given the E(B-V) value quoted in Table 
n the implied N(H)/E(B-V) ratio is consistent with the 
value measured by Bohlin et al. I|1978|l . 



of sight already at low J (and not only at high excitation 
energies) , we present in Fig. |51 a compilation of results 
for H2 levels from J = to 6 for HD 96675, HD 102065, 
HD 10892 7 (Cry et al. lWa . C Oph (Morton [Mil) and C 
Per fSnow irnTTI the Copernicus data have been reanalysed 
using new oscillator strength values as described in Le 
Petit et al. 12004). The first three lines of sight have been 
chosen because the gas probed is known to be far away 
from the star. 

We also include in Fig. results obtained for 
HD 37903, the illuminating star of NGC 2023 in Orion. In 
HST/STIS spectra of that B1.5 V star, Meyer et al. ( 2001 ) 
detect a large amount of highly excited H2, which is con- 
sistent with UV fluorescent excitation of dense molecular 
gas located close to HD 37903. The data shown in Fig. 
El for J < 6 are from ORFEUS observations (Lee et al. 

For the purpose of comparison, column densities have 
been normalised to the total H2 column density appropri- 
ate for each line of sight. It is clear in Fig. that for the 
two lines of sight in which strong absorption from excited 
H2 is detected (from J > 8), HD 34078 and HD 37903, the 
excitation excess is already apparent at 3=2 with respect 
to other, more "standard", lines of sight. Beyond J=3, 
the excess attains two orders of magnitude. In contrast, 
the relative amount in the J=l level (corresponding to T 
w 77 K towards HD 34078) is quite consistent with other 
measurements . 

Apart from the detections presented in this paper, 
several other molecular species have been observed to- 
wards HD 34078 which can be used to constrain a physi- 
cal/chemical model of the gas: 2 (Federman et al. '1994), 
CH, CH+, CN (RoUinde et al. i2003)l and CO (Mac Lachlan 



4.2. Modelling 



To investigate the ability of our model to reproduce ob- 
served column densities we use a recent version of the FDR 
model of Le Bourlot et al. H1993|l as well as the MHD 
shock model of Flower & Pineau des Forets H1998fl . up- 
dated for diffuse clouds and including explicit rotational 
excitation of H2. As described in Le Petit et al. H2002|l . 
the PDR model features a plane-parallel slab of gas and 
dust illuminated on one side by a UV radiation field. 
The radiative transfer as well as the thermal balance are 
calculated (the temperature can also be fixed, if derived 
from observations). The J dependent photodissociation of 
molecules such as H2 and HD as well as their excitation 
in ro-vibrational levels due to radiative excitation, decay 
and coUisional excitation are computed. About 100 chem- 
ical species linked by a network of 1400 chemical reactions 
have been taken into account. The elemental abundances 
and dust size distribution used here are those given in Le 
Petit et al. H2002|l . Note that while the extinction curve 
effectively measured towards HD 34078 is used to describe 
the attenuation of the UV field, the dust size distribution 
adopted is consistent with the average Galactic extinction 
curve. The main properties of the line of sight (extinction 
curve and Ry value) are given in Table Q The density 
and intensity of the incident radiation field in each of the 
two components are adjusted to reproduce the observed 
H2 excitation as well as the amount of other molecules. 
Column densities are computed through integration of the 
local abundances given by the model, assuming a face-on 
geometry. 
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4.2.1. The transluscent cloud 

For the translucent cloud we assume an incident radia- 
tion field equal to that in the solar neighborhood as esti- 
mated by Draine p978|l . An isothermal model is adopted 
at T = 77 K, the kinetic temperature from the two first 
H2 levels. The gas density is well constrained by absorp- 
tion from the first three CI levels. Trying several values for 
iiH, models with densities of 600 cm~^ and below appear 
to underestimate N(CI) whereas those with nn > 800 
cm^"^ overestimate N(CI*). We adopt uh = 700 cm~^ as 
the best compromise. With these parameters, integrating 
abundances to match the observed N(H2, J=0) and N(H2, 
J=l), we find that column densities of all other molecules 
(excepted CH+ and OH discussed below) are well repro- 
duced given the measurement uncertainties and the small 
number of free parameters in the model (cf. Table The 
visual extinction corresponding to this component is Ay 
= 0.8. 

To account for the CH+ abundance we add to the 
PDR component a C-shock with a velocity v = 25 km 
s~^, a pre-shock density of 20 cm^"^ and a magnetic 
field of 7 /iG. The shock contributes only slightly to 
the excitation of H2 for levels J = 5 to J = 7 (Fig. [T]). 
The thickness of the shocked layer corresponds to Ay ~ 
0.1 and gives a negligible contribution for other species, 
OH excepted. In this latter case, the contribution of the 
MHD shock is too large by a factor of about 7 (Table 01 . 
The maximum temperature of neutral species reached in 
the shock is 4000 K. A single shock has been used here; 
vortices or a collection of several weaker shocks could 
have been considered as well to represent the dissipation 
of kinetic energy along the whole line of sight (cf. Gredel 
et al. I2UU2|I . which might help to reduce the discrepancy 
between the observed and model N(OH) values. However, 
since the C-shock required to account for N(CH+) 
appears to play no important role for the H2 excitation, 
no further attempt to fit both the observed OH and CH+ 
values was made (OH observations are relatively scarce 
due to the wavelength range of its absorption lines and 
thus little is known about correlations between OH and 
CH+). 



4.2.2. The hot PDR around HD 34078 

The contribution of material surrounding HD 34078 is 
strongly suggested not only by the anomalous excitation 
diagram but also by the IRAS detection of a far infrared 
excess associated with this star by van Burenet al. l|TM5j) . 
According to these authors, such an emission is indica- 
tive of a bow shock located around the star, at the inter- 
face between the stellar wind and the ambient interstellar 
medium, where the material is strongly compressed. In 
this location, dust grains are exposed to an intense UV 
field and emit FIR radiation, accounting for the excess 
observed in IRAS data. We shall assume in the following 
that a scenario like the one proposed by Mac Low et al. 



H1991II to describe cometary compact HII regions also ap- 
plies in our case and that H2 is present in a thin layer 
within the bow shock (cf their Fig. 2). This is not obvious 
here since the ambient medium is not necessarily dense 
enough to contain H2 molecules (this point is discussed 
in Sect. 4.2.3). Below, we consider this bow shock as a 
PDR and show that excitation of the H2 molecules by UV 
photons can explain the observed excitation. 

Since we have no observational constraints on the tem- 
perature for this component, we now consider models in 
which the thermal balance equation is solved. The den- 
sity n and the intensity of the incident radiation field x 
expressed in Draine's units are varied so as to reproduce 
the H2 excitation diagram. We find that adopting tih= 
10^ cm~'^, X = 10'' and integrating abundances up to 
Ay — 1.0, we get a remarkably good fit given the lim- 
ited number of free parameters, as illustrated in Figs. [SI 
and|3 H2 column densities for levels J = 3 and above are 
correctly reproduced, indicating that the hot PDR natu- 
rally accounts for both the overpopulation measured for 
J > 3 and the high energy tail (beyond J = 7) not seen 
in most lines of sight (let us recall that N values are well 
constrained by observations both for J = 3 and J > 8 since 
in these two limits W no longer depends on b). Then, the 
"standard" H2 J = 2 to 7 population as shown in Fig. El 
is masked on this particular line of sight by the hot PDR 
component. Most of the observed atomic hydrogen comes 
from this component. Because of the high radiation field, 
column densities of molecules like CO, CH, OH etc are 
negligible as compared to those of the translucent compo- 
nent due to large photodissociation rates (cf. table 0J. 

By varying n and x around the best fit values, we find 
that the density is well constrained by the model: with 
X — W^, the amount of high excitation H2 can be re- 
produced only with uh between 5000 and 5 10** cm^"^. 
On the other hand, the intensity of the radiation field is 
poorly constrained. For a density uh = 10^ cm~'^ it is 
possible to reproduce the H2 excitation diagram without 
exceeding the observed N(H I) and N(H2) with x between 
500 and 10**. We thus find that a two-component model 
(complemented by a C-shock to account for CH+ produc- 
tion within the translucent cloud) successfully reproduces 
the main characteristics of the gas along the line of sight 
towards HD 34078. 

A comparison of the excitation diagrams obtained for 
HD 34078 a nd HD 37903 (Fig. 6, Fig. 7 and Fig. 2 from 
Meyer et al. I2()()l|l indicates that they are relatively simi- 
lar. Moreover, it is noteworthy that the amounts of H2 and 
dust are comparable on these two lines of sight: N(H2) = 
6.4 10^0 and 7.8 10^° cm^^ and E(B-V) = 0.52 and 0.35 
for HD 34078 and HD 37903 respectively (the spectral 
types are also close to each other: 09.5V and B1.5V). It 
is then not surprising that when modelling the high en- 
ergy part of the HD 37903 diagram, Meyer et al. H2001|l 
derive values for the radiation field and gas density consis- 
tent with our estimate for the hot PDR component around 
HD 34078. In our picture, the markedly different proper 
velocities of these two stars (only HD 34078 is a runaway) 
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Fig. 7. Observed H2 excitation diagram and H2 column 
densities for each of the three components in the model. 

could just correspond to distinct mechanisms responsible 
for the presence of high density gas close to each star. 

4.2.3. Discussion 

As already noted, we adopted in the above model the true 
HD 34078 extinction curve but a dust size distribution 
corresponding to the average Galactic extinction curve. 
This lack of consistency could have significant effects on 
model results. To estimate the uncertainties related to the 
(unknown) dust size distribution, we ran additional mod- 
els using the average Galactic extinction curve (including 
the appropriate Ry value) instead of the HD 34078 one 
and imposing similar values as those given in Table 01 for 
the H2 column density of the translucent cloud and hot 
PDR components. We find that the results of both models 
remain within the observational uncertainties; in the H2 
excitation diagram, the maximum difference in N(H2, v, J) 
is less than 25% and similarly small relative variations are 
seen for the abundances of species like CI, CH, CN (the 
variation attains 50% for CO). We then conclude that our 
model results are not strongly dependent on the adopted 
dust size distribution and that the good agreement with 
observations is fully significant. 

In order to better assess the reliability of the hot PDR 
component, one can estimate from dynamical arguments 
the distance of the bow shock to the star and infer the 
corresponding radiation field intensity. Using Eq.(2) from 
van Buren et al. (1990), adopting an average density of 
the ambient ISM of Ua ~ 1 cm"'^ and a star velocity of 
100 km s^^, we get a distance of dg ~ 0.5 pc. Considering 
the observed flux of HD 34078 at A ~ 1100 A (10"" 
erg.cm~^.s^^.A~^) and correcting for the dust attenua- 
tion implied by E(B-V) — 0.5 and the observed extinction 
curve (thoo — 6.3), we find that the distance drad at which 
the intensity of the radiation field of the star is equal to 
10^ times the ISRF is 0.2 pc, approximately the estimated 
distance of the shock to HD 34078. This good agreement 



gives coherence to our model. An additional test would 
consist in estimating the density of the material located 
in the bow shock, but unfortunately this is strongly depen- 
dent on the magnetic field and the shock velocity (both 
unknown) . 

Although the spectral resolution is limited, one can 
also use the absorbing gas kinematics to test our model. 
In the latter, all H2 absorptions from levels above the 3=2 
one arise from the same component, the hot PDR. These 
lines should therefore be characterized by a common ve- 
locity. To search for a possible Doppler shift among J > 3 
lines we selected spectral regions with good S /N ratio and 
displaying J = 3, 4, 5 lines with closeby J = 9 absorptions. 
Similarly, around 1059 A, a group of 7 H2 lines from v ~ 
0, J = 4, 6, 7, 8 and v = 1, J = 3, 4 is present. In all 
these cases, we find that the r.m.s. scatter in wavelength 
shifts is no larger than 0.005 A, which implies a velocity 
difference smaller than 1.5 km s~^ (or about 5 km s~^ at 
3cr). 

Conversely, since the translucent and hot PDR compo- 
nents are unrelated in our model, their velocity distribu- 
tions should differ both in the average and width, unless 
they happen to coincide by chance. J = and 1 absorp- 
tions are too broad to provide any valuable velocity esti- 
mate. We can use instead HD lines which trace the same 
component. Near the 1031 and 1054 HD J = absorptions, 
we find a few narrow H2 absorptions (from the hot PDR 
component) and by comparing measured wavelengths to 
the laboratory values, we find no significant shift between 
the two components to a similar limit as above (5 km 
at 3it). The test is however not so constraining because 
the number of measurements is limited and the accuracy 
of the latter is difficult to assess because i) the continuum 
is uncertain and ii) unrecognized blends may be present. 
Regarding velocity dispersions, one would expect a pri- 
ori a larger value for disturbed gas near HD 34078 (i.e. 
the hot PDR component) than for the foreground translu- 
cent cloud. High resolution observations of CH absorption 
can be used to probe the latter, which indicates b w 3 
km s-i (Allen EMI RoUinde et al. I^HMjl . As discussed 
in Sect. 2.2, there is no clear evidence for a significantly 
larger dispersion in high excitation gas. Thus, tests involv- 
ing kinematics do not provide additional support for our 
two-component model. 

The similarity of the velocity and velocity dispersions 
of the two components suggests that perhaps the latter 
are not unrelated. The suggestion made by Herbig H1958|l 
is interesting in this regard. Based on a comparison of Ha 
and visible continuum images, he proposed that HD 34078 
has recently encountered new cloud material. In such a 
scenario, one can imagine that the density of the ambient 
medium, Ua, is about 10-20 cm"'^ instead of 1 cm^^. A first 
consequence involves the origin of the excited H2. With 
Ua — 1 cm~'^, H2 molecules would not preexist and would 
need to be formed in the bow shock; we estimate that at n 
~ 10^ cm^"^, 10*^ yrs are required while it takes no longer 
than lO'' yrs (a few 0.2 pc/100 km s^^) for the compressed 
material to drift from the head of the shock to the sightline 
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(after lO** yrs, the fraction of molecular H2 is less than 10% 
of its steady-state value). In contrast, if Ua — 20 cm^^, 
H2 can be supposed to be already in place in the ambi- 
ent medium. Moreover, applying mass conservation for the 
swept gas leads to a predicted H2 column density through 
the hot PDR which is consistent with the observed one. 
Indeed, one gets the relation TrR^.Ua = 2-K.R.A.R.npDR 
or Ua-R = 2.AR.npDR — 2Np£iR where R is the radius 
of the cylinder containing the swept material (which we 
assume to be comparable to ds), Ai? the thickness of the 
PDR, npDR and Nrdr the density and column density 
throughout the compressed layer. We find that Nrdr is of 
the order of 10^^ cm"^, comparable to the column density 
estimate for the hot PDR drawn from our best fit model 
(Table Finally, in a denser ambient medium the shock 
will occur closer to the star {dg goes as n~'~'-^) and with 
Ua — 20 cm~'^, ds ~ 0.1 pc (instead of 0.5 pc for 1 cm^'^), 
quite comparable to drad- Thus, a relatively high density 
for the ambient medium around HD 34078, as suggested 
by Herbig f 1958)1. appears to be well consistent with the 
modelling of the II2 excitation. 

One question that remains open is the relation between 
the translucent component containing the bulk of the H2 
(J = 0, 1) and the hot PDR material. The low temperature 
(77 K) of the II2 gas, the low velocity dispersion as esti- 
mated from CH and the presence of many molecules that 
might be photodissociated if located too close to HD 34078 
suggest that the translucent cloud is not associated to the 
star. In contrast, the similarity of velocities and velocity 
dispersions favors a picture in which the two components 
have a common origin. Below (Sect. 6), we propose a few 
observations which should help to clarify the relation be- 
tween HD 34078 and the bulk of the molecular gas. 

5. Absorption line variations and small scale 
structure towards HD 34078 

The original purpose of our project was to probe the 
structure of quiescent molecular gas. The discovery of 
high excitation H2, making the line of sight towards 
HD 34078 somewhat peculiar, could have been embar- 
rassing. However, the modelling work described above in- 
dicates that absorption from the J = and 1 H2 levels 
(which comprise by far most of the gas) traces material 
unperturbed by HD 34078. In the following, we first search 
for variations in the damped line profiles from these two 
levels, which can be used to investigate in a direct and 
efficient way column density variations, regardless of any 
assumption on the velocity distribution. The five spectra 
available span a time interval of 2.7 yr, which corresponds 
to linear scales of about 50 AU in the foreground cloud. 

In our model, lines from all other H2 levels originate es- 
sentially from gas surrounding HD 34078. Since the com- 
pressed material is continuously flowing across the line 
of sight as the star moves, the time behavior of these 
lines should tell us whether this flow is regular or chaotic. 
Further, if there exist density fluctuations in this gas or 
changes in the UV field, this should also result in varia- 



tions of absorption lines arising from the hot PDR com- 
ponent. In order to select those H2 levels that display the 
highest sensitivity to changes in n or x, we use the model 
presented in Sect. 4. We vary the density, n, and radia- 
tion field, X, (separately) by a factor of 2, and find that 
the sensitivity to changes in x is small for all levels (30 % 
at maximum for J = 4 - 5), while the response to similar 
variations in n is much larger with a peak at J = 6. N(J = 
6) varies by a factor of 3.4, and beyond J = 8 the induced 
variation in N(v, J) is around 50%, relatively independent 
of the level. Unfortunately, absorptions from levels J = 5, 
6, 7, 8 are heavily saturated and the observed profiles or 
W values are weakly dependent on N(J) (these lines can 
be used instead to assess the stability of the instrument). 
We shall then examine the good S /N lines with intermedi- 
ate opacity (r w 1 — 3) arising from high excitation levels 
(J>9). ' 

In Boisse et al. (|l2"0T)!i)l . we presented a report based on 
the first three spectra, among which no significant changes 
were noted, except a possible variation of lines from highly 
excited H2. 

5.1. Sensitivity of absorption lines to column density 
variations 

The sensitivity of an absorption line profile to small col- 
umn density variations can be estimated from the ratio 
Q = ((5///c)/((5N/N), where / is the o6serwerf profile and Ic 
the continuum. Considering ideal Voigt profiles observed 
at a given finite resolution, Rollinde l|2()()3|l has shown that 
a good sensitivity is achieved in two cases: i) fully resolved 
lines (e.g. damped H2 absorption from low J levels) and 
ii) marginally thin lines (central opacity, tq ~ 1). 
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Fig. 8. Profile variation for the H2 Lyman (4-0) R(0) 
1049.37 A transition. Top panel: The red line corresponds 
to log(N) = 20.6. The dashed profile is obtained by vary- 
ing N by -t- 10%. Lower panel : difference between the 
two profiles. In both panels, regions noted (a and b) are 
subject to little variation which implies the presence of a 
maximum (M) in between. 
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Consider an unblended and fully resolved damped line. 
If the column density increases the profile broadens, which 
results in a decrease of I/Ic at any given wavelength, ex- 
cept in the line core and outside the wings where I = Ic 
(Fig.lS)). It can be shown that the Q parameter reaches its 
maximum, Qmax = 0.357, on each wing of the profile at 
wavelengths where the optical depth tq equals 1 and I /Ic 
= 0.37 (RoUindeHnnSI)- For a change of N by 10%, the 
expected variation in I/Ic is 0.035 (well above the noise 
level in our FUSE spectra). The corresponding broaden- 
ing of the profile at tq = 1 on each wing is typically 0.05 
A (a value comparable to the resolution). 

A damped profile corresponds to 

"6.726 x 10-1^1 N/A^7 
_ J (Ao/A-1)^ 

where c is the light velocity (in km s^^), f the oscillator 
strength, 7 the radiative rate (in s"^), and where Ao and 
N are expressed in A and cm~^ respectively. Then, if vari- 
ations are observed for one specific line, one can predict 
exactly the corresponding variations expected for other 
damped lines from the same level. This gives a powerful 
method for checking the reality of any variation and rule 
out artefacts of instrumental origin (e.g. due to distortions 
in the wavelength calibration). 

In practice, the low J II2 transitions in FUSE spectra 
appear in "systems", where one J = 0, two J = 1 and 
one 3 = 2 lines are strongly blended (e.g. for the (4-0), 
(3-0) and (2-0) Lyman bands at A « 1050 A, 1063 A, 
1078 A respectively). To compute the expected variation 
of the whole system profile, we assume that the excitation 
of the gas is uniform i.e. that the relative variations of N(J 
= 0), N(J = 1) and N(J = 2) are identical. Comparison 
of synthetic profiles corresponding to values of N(J = 0, 
1 and 2) similar to those given by profile fitting of the 
FUSE spectra indicates that the blue wing of the J = 
line and the red wing of the 3 = 2 line in each system (the 
two extreme transitions of the blend) are relatively little 
affected by the J = 1 transitions (Fig. inj. For the 1050 
A H2 system, changes in N(J = 1) alone result in profile 
variations that would not be easy to detect because they 
appear in the far wings of the blend; as a consequence, no 
useful constraint can be obtained for N(J = 1) from that 
system. 

5.2. Comparison of H2 clamped profiles obtained 
between Januray 2000 and October 2002 

To search for variations among the observed spectra one 
could fit each system profile (or fit them simultaneously) 
and then compare the N(J=0, 1, 2) values obtained for 
each epoch. However, the accuracy of individual N values 
would then be strongly limited by our knowledge of the 
stellar continuum, which is poor given the large width of 
the II2 absorptions and high number of additional absorp- 
tion lines present. A great advantage of reobserving the 
same star is that variations can be searched for through a 
direct comparison of the spectra themselves. 
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Fig. 9. Same as in Fig. |Slbut for the whole blend of J = 
0, 1, 2 lines ((4-0) Lyman band). The assumed variation 
is d N/N = 10 % for each of these three levels. 



Before performing this comparison, we need to adjust 
the wavelength and flux scales of the two spectra consid- 
ered. These two tasks are performed separately for each 
H2 system, on an extract of the spectra centered on the 
absorption and covering about 15 A. Using narrow absorp- 
tion lines located on both sides of the broad absorptions, 
the wavelength shift to be applied is measured. Similarly, 
a multiplicative factor bringing the adjacent continua of 
the two spectra to the same level is determined. In fact 
the "adjacent continuum" is not easy to define because 
the various systems overlap, and given our N(H2, J) val- 
ues for J = 0, 1,2 the entire spectral domain of interest is 
affected by H2 absorption. Thus, to perform the relative 
flux normalisation, we use regions where H2 absorption is 
minimum, about 7 A away from the center of each trough 
(e.g. near 1043, 1057, 1071 A). It is to be noted that 
such a procedure is approximate only, because a varia- 
tion in N(Il2) will slightly affect the continuum level, and 
imposing the same continuum value in the two spectra 
to be compared will result in an apparently lower vari- 
ation. However, in our case, the absorption in between 
two H2 systems is no larger than about 10%. By apply- 
ing the inter-normalisation procedure described above to 
synthetic spectra with N(H2, J) differing by 10%, we find 
that the apparent reduction in the N variation is not re- 
ally significant, given our accuracy. In the following, we 
therefore ignore this difhculty 

These two corrections (involving two values only for 
a given H2 system and pair of spectra) account i) for un- 
known changes in the zero-point of the FUSE spectra from 
one epoch to another (and possibly, for inaccuracies in 
the FUSE wavelength calibration) and ii) for stellar flux 
variations (and/or variations in the flux calibration). The 
second correction assumes that no variation in the spec- 
tral shape of the continuum has occurred, which appears 
consistent with the data (except blueward of the OVI dou- 
blet where impressive variations of the absorption due to 
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the stellar wind are seen) . As shown in Fig. IIUI the super- 
position is excellent once the two corrections are applied, 
which indicates a very good instrumental stability. 

Variations of the H2 systems among all spectra have 
been examined (see Boisse et al. 120011 for a comparison 
of the first and third spectra). Here, we present results 
involving the first and fifth spectra, which provide the 
largest timebase. The comparison has been performed for 
each of the three 1050, 1063 and 1078 A H2 systems. In 
Fig. ^1 we can see that the difference near the red or blue 
end of the through where the largest variations are ex- 
pected is no larger than a few percent of the continuum 
(windows of width 0.3 A are considered) and that no 
systematic trend towards an increase or decrease of N(H2) 
is noticeable (the curves shown in the lower panel are for 
N(H2) changes by ± 10 %). For systems other than those 
around 1050, 1063 and 1078 A, the S/N ratio is too low or 
the presence of additional absorption strongly alters the 
H2 profile in the spectral regions where the largest varia- 
tions are expected for I / Ic ; in particular, we exclude the 
1037 A system which is affected by variable OVI absorp- 
tion from the stellar wind. 

Very recently (September 2003), we obtained a sixth 
FUSE spectrum which brings the time interval and scale 
covered to 3.7 yr and 63 AU respectively. Although these 
data have not been fully analysed yet, we could verify 
that the profiles of the three H2 systems considered above 
are again completely consistent with those of the earlier 
epochs. 



1045 1050 1055 




-0.05 -■' 



1045 1050 1055 

>.(A) 

Fig. 10. Comparison of the January 2000 and October 
2002 HD 34078 FUSE spectra in the region comprising 
the 1050 A H2 system. Top panel: Superposition of the 
normalised spectra (black line: Jan. 2000; red line: Oct. 
2002), after applying the corrections discussed in the text. 
Lower panel : The difference between the two observed 
normalised profiles. The solid lines indicate the variation 
expected of the H2 absorption profile for a change (5N/N 
= ± 10% of the column density for each of the J = 0, 1 
and 2 levels. 



We thus detect no appreciable change in the damped 
H2 profiles. To get an upper limit on column density vari- 
ations, we shall use a Maximum Likelihood Estimator (as 
described in Numerical Recipes, Sect. 15.1). The data cor- 
respond to the flux variations along the profile of the three 
systems mentioned above, while column density variations 
are the parameters to be determined. For simplicity, we 
shall consider that H2 excitation among the J = 0, 1 and 
2 levels is uniform, i.e. constant in time (clearly, local 
changes in the excitation would constitute a second or- 
der effect). In this assumption, the observed profiles that 
are blends of absorption lines arising from these three lev- 
els depend only on one single parameter, the total (or 
relative) N variation. 

The likelihood of SN given a flux variation SI is a func- 
tion of the wavelength, written as P{SN\5I, A). This prob- 
ability distribution is computed from random realisations 
of synthetic spectra with different column densities (us- 
ing our fit to the data as a reference) to which we add 
white noise with an r.m.s. scatter similar to that mea- 
sured in the observed spectrum. Then, we get the com- 
bined likelihood function P{SN\5I) for the set of three H2 
systems quoted above by computing the normalised prod- 
uct of the three individual distributions. The data imply 
SN/N = 0.02 ± 0.01, which corresponds to 3cr lower and 
upper bounds, -1% < SN/N < 5%. 

Observations of the CH and CH"*" blue lines have 
suggested possible correlated variations of N(CII) and 
N(CH+) (Rolhnde et al. I2003|l on the scale of about one 
year. One may then wonder in case these fiuctuations are 
real (additional observations are in progress to check this) , 
whether they are accompanied by similar N(H2) varia- 
tions. It turns out that, for two CH measurements which 
display the largest variation in W(CH A4300) {+ 6.3% be- 
tween January 2000 and February 2001 corresponding to 
-I- 9.5% in N(CH)), there exist nearly simultaneous FUSE 
observations (to within a few days). We then examined 
carefully the 1050 A H2 system (which provides the tight- 
est constraints) and searched for any hint of an increase 
in N(H2) between the two epochs. No such change can 
be seen and by comparing synthetic profiles obtained by 
slightly varying N(H2, J = 0, 1, 2) about the best fit val- 
ues we find that an increase of N(H2) as large as 9.5% is 
clearly ruled out by our data. 

Let us conclude this section by a few comments on 
some systematic effects that might affect the analysis de- 
scribed above. For instance, slight distortions in the wave- 
length scale could, if large enough, mimic the kind of sig- 
nature we are searching for (a shift of the blue or red end 
of the H2 troughs). Local fiuctuations in the continuum 
level due to varying stellar absorption lines could also in- 
duce apparent H2 absorption changes. However, among 
the three systems considered these effects are found to be 
negligible, as indicated by the good agreement between 
the wavelength shifts derived from various narrow lines 
adjacent to each system. If future observations were to 
suggest real variations, we note that the spectral shape of 
the expected variation could be used to recognize these ef- 
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fects; moreover, data from other detectors (such as LiF2B, 
which covers in part the same range as LiFlA) could be 
used also (although of lower S /N ratio) to rule out instru- 
mental effects. 

5.3. Lines from highly excited H2 

We now search for variations among the low optical depth 
lines arising from high H2 levels and tracing (in our model) 
the hot PDR component around HD 34078. In order to 
avoid false detection of variations due to instrumental 
changes (e.g. spectral resolution) or to variations in stellar 
features blended with the H2 line considered, we require 
that at least two lines from the same level display the same 
behavior. We then focus on those levels from which several 
lines with a high degree of significance and intermediate 
opacity values (tq ~ 1 - 5) are detected (e.g. J = 9, (v = 1, 
J = 1, 2, 3)). Occasionally, we observe apparent variations 
e.g. for the V = 1 lines at 1052.6 A (J = 1) and 1056.8 A 
(J =3) for which a decrease in W by about 20 % between 
epoch 1 and 2 is suggested by the data. However, we can- 
not find any H2 level for which unambiguous evidence for 
a column density variation can be obtained. Often, the 
shape of the continuum adjacent to absorption lines is not 
stable (probably because of changes in the profile of un- 
derlying stellar absorption) which greatly complicates the 
comparison of spectra and the estimate of uncertainties 
on W or N variations. We thus conclude that the FUSE 
data exclude the existence of large variations in excited H2 
lines (larger than about 30%) . Clearly, HST spectra would 
be very helpful to better constrain the temporal behavior 
of highly excited H2 absorption thanks to both a higher 
resolution and S/N ratio. 

5.4. HI 

The Ly/3 line offers an interesting opportunity to search for 
N(HI) variations. As already mentioned, the N(HI) value 
obtained from modelling the H2 - Ly/3 blend is larger than 
the estimate drawn from lUE observations (Mc Lachlan & 
Nandv [T984|l . In Fig. ^2 we show the spectral region used 
to determine N(HI) from the Ly/3 - H2 blend together with 
synthetic profiles computed with the 1979 and 2000 N(HI) 
values. 

To assess whether the two determinations are really 
inconsistent, we have reanalysed the Lya data. In the 
spectrum that we retrieved from the lUE database there 
is some weak residual flux at the bottom of the damped 
Lya profile (while, given the spectral resolution, the pro- 
file should reach the zero level in the absence of scattered 
light). Correcting for this offset and fitting the Lya pro- 
file again, we get a value - log(N(HI)) = 21.30 - slightly 
larger than that given by Mc Lachlan & Nandy fl984): 
log(N(HI)) = 21.24. This revised N(HI) remains signifi- 
cantly lower than the January 2000 value. Although it is 
difficult to reach a firm conclusion regarding the reality 
of a change in N(HI) because we could not compare the 



same feature from both epochs and because the Ly/3 line 
is strongly blended with H2 absorption, our analysis sug- 
gests that the HI column density towards HD 34078 has 
increased by 38% between 1979 and 2000 (this time in- 
terval corresponds to a transverse separation of 300 AU if 
the HI is assumed to lie at d = 400pc) . The corresponding 
average increase rate for N(HI) would be 1.8%/yr, which 
turns out to be comparable to the rate implied by CH 
observations (1.7%/yr; RoUinde et al. lM)^ . 




1024 1026 
Wavelength (A) 



Fig. 11. The bottom of the H2 - Ly/3 trough near 1025 
A. The normalized January 2000 (thick black line) and 
October 2002 (thin red line) are shown together with syn- 
thetic profiles obtained assuming N(HI) = (long-dashed 
red line), log(N(HI)) = 21.30 (the value inferred from 1979 
lUE Lya data; short-dashed black line) and log(N(HI)) = 
21.44 (our best value; dotted blue line). Note how the lat- 
ter value improves the fit on both edges of the trough. The 
line at 1023.4 A is from H2 J = 4. The N(H2) values for 
J = 0, 1, 2, 4 used to compute profiles are those given in 
Table 1^21 The emission feature cut out near 1025.5 A on 
both spectra is terrestrial Ly/3. 



On smaller scales, the ubiquitous structure in the 
atomic phase (Heiles |1997 ) should be apparent if our 
sensitivity to variations in the Ly/3 profile is good enough 
(most of the HI is likely unassociated with HD 34078). 
Modelling of the profile indicates that the region near 
1026.9 A (rest wavelength), at the red end of the trough, 
is best suited to search for N(HI) variation. At longer 
wavelengths, where larger variations would be expected, 
strong lines are present as well as OVIA1031 variable ab- 
sorption from the O star wind. We carefully compare the 
five profiles in this region (after adjusting the wavelength 
and flux scales as explained in Sect 4.2 above) and flnd 
no evidence for any variation larger than 10% (the exact 
upper limit is difficult to define but we estimate that 
a 10% variation would unambiguously be seen; 10% is 
therefore equivalent to a 3 cr limit). 
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5.5. Discussion 

The main goal of our work is to study the AU scale struc- 
ture in the spatial distribution of H2 through the tempo- 
ral behavior of its absorption lines. Concerning the bulk 
of molecular gas (i.e. J = and 1 levels) which is be- 
lieved to lie far from HD 34078, no profile variation is 
detected, which implies column density fluctuations lower 
than 5% over scales ranging from 5 to 50 AU. We conclude 
that translucent clouds display no " ubiquitous" small scale 
structure in molecular hydrogen such as that seen for 
H2CO, OH or in the atomic phase through HI 21 cm ab- 
sorption (Heiles .1997)) or Nal optical lines. Our sampling 
remains limited, but since we observed five distinct posi- 
tions, from 5 to 50 AU apart, one can reliably conclude 
that a line of sight with N(H2) ~ 3. x 10^" cm^^, does not 
in general intersect clumps of this length scale with N(H2) 
larger than 1.5 x 10^^ cm~^. 

The observed absence of ubiquitous N(H2) fluctuations 
in translucent clouds is an important result in that it tells 
us that uniform models can legitimately be used to de- 
scribe chemical and physical processes within interstellar 
gas (which we did above in Sect. 4 !). Structure might be 
present at larger scales and it is noteworthy that direct CO 
emission imaging using millimeter wave interferometers 
can now probe structure at scales as small as 500 AU, i.e. 
barely one order of magnitude larger than those explored 
here. A few projects of this type are in progress; prelimi- 
nary results indicate some velocity structure but modest 
column density variations (Falgarone, private communica- 
tion). 

CH observations have indicated a 20% increase of 
N(CH) over the last 12 years (Rolhnde et al. lMI^ . Since 
CH and H2 appear to be well correlated from one line 
of sight to another, one can search for any such trend 
in the FUSE data. The two extreme FUSE observations 
are separated by 2.7 yr and show no hint of such a varia- 
tion. Assuming that the rate of increase is constant (about 
1.7%/yr for N(CH); in fact the shape of this variation is 
not well constrained: cf RoUinde et al. I2()()3|l . an increase 
of N(H2) by 4.5 % would be expected if a hnear correla- 
tion still holds at small scales. This is just below what we 
can detect. Therefore, the present data do not allow us to 
check whether the long-term evolution of CH and H2 is 
similar. 

Regarding highly excited H2, the time behavior of its 
absorption lines is related, in the frame of the model 
sketched in Sect. 4, to the stability of the gas flow in the 
shocked gas layer, at the interface of the stellar wind and 
the ambient medium in which HD 34078 is moving. If the 
latter is not homogeneous in the 1 - 100 AU scale range, 
this will induce line variations resulting from fluctuations 
of the amount of gas swept up by the bow shock and 
then compressed in the cylindrical layer around HD 34078. 
Variations due to instabilities in the interface could also 
occur, even if the ambient medium is homogeneous. The 
lack of pronounced variations indicates that the medium 



through which HD 34078 is moving is smooth and, further, 
that marked instabilities do not develop at the interface. 

6. Summary and Prospects 

The analysis of five HD 34078 FUSE spectra taken be- 
tween January 2000 and October 2002 leads to the follow- 
ing conlusions: 

- the unexpectedly large amount of higly excited H2 
is attributed to the presence of molecular gas in a dense 
layer surrounding HD 34078 and excited by fiuorescence. 
Detailed modelling of this region provides a very good fit 
to the excitation diagram, with few free parameters, 

- the bulk of the molecular gas (H2 (J = 0, 1), CH, 
CH+, CO etc) lies in a foreground translucent cloud, pre- 
sumably unrelated to HD 34078. Physical conditions (tem- 
perature, density, etc) are similar to those seen in diffuse 
molecular gas in general, 

- the stability of damped H2 absorption lines from J = 
0, 1 and 2 levels indicate that H2 molecules are smoothly 
distributed within the foreground translucent cloud with 
column density fluctuations lower than 5% over scales 
ranging from 5 to 50 AU. 

- H2 lines from highly excited levels also turn out to be 
stable. Again, this points towards a smooth distribution 
of H2 molecules in the diffuse ambient medium through 
which HD 34078 is presently moving and indicates that 
no instability is developing in the shocked gas flow around 
HD 34078. 

Let us conclude with a few prospects concerning the 
study of this very interesting line of sight. As discussed 
above, our analysis of small scale structure assumes that 
the translucent cloud lies far away in front of HD 34078 
and thus heavily relies on a good understanding of the lo- 
cation of the gas responsible for the various absorptions 
observed. It is then important to find additional support 
for the two-component model proposed. We first plan to 
map the CO emission in the close environment of the 
HD 34078 line of sight. If our assumption is correct, there 
should be no relation between the morphology of the ob- 
served CO emission and the position of the O star or the 
direction of its motion (essentially northward) . While com- 
paring the radial velocity (and velocity dispersion) of the 
absorptions arising from each of the two main components, 
the limited FUSE resolution prevented us to reach definite 
conclusions. HST/STIS observations would be invaluable 
in this respect. Thanks to a higher resolution and S/N, 
HST spectra would also allow us to much better charac- 
terise both the kinematical and physical properties of the 
translucent and hot PDR components, through detection 
of H2 lines from higher excitation levels (cf the study by 
Meyer et al. FSDOl,) as wefl as from CI or CO. Finally, the 
Lya profile would give an accurate N(HI) determination 
and then allow us to confirm the reality of the long-term 
N(HI) variation. 
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Table 2. Measured column densities for various H2 levels (best estimate and 3 a lower and upper bounds) and 
predictions from our model (see text). The number of lines effectively used in the determination of N(H2, v, J) is given 
as Ml 
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Table 3. Equivalent widths and column densities derived from the OH transitions observed towards HD 34078. 
Wavelengths and oscillator strengths are taken from Felenbok and Roueff (1996). 
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Table 4. Observed and computed column densities in cm 
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1.2-10" 


1.7-10^^ 


1.7-10^^ 




6.4-10^" 


6.0-10^° 


6.9-10^" 


6.5-10^° 


5.1-10" 


2.7-10" 


6.8-10^° 


Av 


1.8 






0.8 


0.1 


1.0 


1.9 


HD(2) 




10" 


10" 


8.5-10" 


— 


5.4-10" 


8.5-10" 




9.5-10^^ 


9.3-10" 


1.1-10" 


5.0-10" 


4.1-10" 


7.4-10^ 


5.4-10" 




6.5-10^^ 


5.8-10" 


7.1-10" 


1.9-10" 


6.0-10" 


5.0-10" 


6.0-10" 




5.8-10" 




— 


2.3-10" 


2.9-10" 


1.4-10^ 


2.3-10" 




2.2-10^^ 


1.7-10^^ 


2.7-10^^ 


7.2-10" 




1.3-10^ 


7.2-10" 


OH(2) 


3.5-10" 


1.4-10" 


5.6-10" 


1.4-10" 


2.5-10" 


5.9-10" 


2.6-10" 




3.3-10^2 






1.6-10" 


5.8-10" 


4.8-10* 


2.2-10" 




5.7-10" 


4.6-10" 


7.2-10" 


1.1-10" 


3.0-10" 


1.1-10" 


1.1-10" 


CI(2) 


9.4-10" 


3.6-10" 


1.7-10" 


3.5-10" 




2.4-10" 


3.5-10" 


cr(2) 


5.8-10" 


1.6-10" 


5.8-10" 


5.7-10" 




6.9-10" 


5.7-10" 


QJ**(2) 


2.2-10" 


l.l-lO" 


4.0-10" 


2.8-10" 




i.o-io" 


2.8-10" 



Dash means an unknown value. 

Last column (Total) corresponds to the sum of the components in the model. 
References for the observations are: 



(1) This work and Mc Lachlan & Nandy (1984). 

(2) This work. 

(3) RoUinde et al. (2003). 

(4) Federman ct al. (1994), RoUinde et al. (2003) 

(5) Oka et al. (2003). 

(6) Mc Lachlan & Nandy (1984) 



